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ABSTRACT 

The luminosities of SGRB host galaxies are anticorrelated with both the isotropic equivalent gamma 
ray energy and the gamma ray luminosity of the explosions. Observational selection effects only 
strengthen the significance of this correlation. The correlation may indicate that there are two phys- 
ically distinct groups of SGRBs. If so, it requires that the more luminous class of explosions be 
associated with the younger class of progenitors. Alternatively, it could be due to a continuous dis- 
tribution of burst and host properties. As one possible explanation, we find that the effect of binary 
neutron star masses on inspiral time and energy reservoir produces a correlation of the appropriate 
sign, but does not automatically reproduce the correlation slope or the full range of SGRB energy 
scales. Any future model of SGRB progenitors needs to reproduce this correlation. 

Subject headings: gamma-rays: bursts 



1. INTRODUCTION 

The host galaxies of short gamma ray bursts (SGRBs) 
exhibit a wide range of physical properties. SGRBs have 
been found not only in star-forming hosts, but also in 
elliptical galaxies with strong upper limits on their star 
formation rate (hereafter "quiescent hosts") (Prochaska 
et al 2006; Berger 2008). In contrast, long gamma ray 
bursts (LGRBs) are found exclusively in star forming 
galaxies (e.g., Le Floc'h et al 2003, Fruchter et al 2008, 
Savaglio et al 2009). Several lines of evidence now link 
LGRBs with the deaths of massive stars (Hjorth et al 
2003, Stanek et al 2003, Fruchter et al 2006, Raskin et al 
2008). The origin of the SGRBs is not firmly established, 
but suspicion centers on merging neutron star binaries, 
or other phenomena primarily associated with neutron 
stars (Belczynski et al 2006; Chapman, Priddey, & Tan- 
vir 2008). Hosts offer one of the best tools for unveiling 
the nature of the GRB progenitors, because the prop- 
erties of the gamma ray emission and afterglow are de- 
coupled from many important details of the progenitor 
objects under the fireball model of GRBs (e.g. Paczyhski 
& Rhoads 1993, Katz 1994, Piran 2005). 

The situation is similar to supernovae, which divide 
into two groups, the type la and the core collapse events 
(which include types II, lb, and Ic). Core collapse events 
are found only in star forming galaxies, and are asso- 
ciated with the deaths of stars having initial masses 
> 8Mq. In contrast, type las are found in both qui- 
escent and star forming hosts, and are thought to be 
powered by the nuclear explosion of a white dwarf whose 
mass is pushed above the Chandrasekhar limit by accre- 
tion. However, despite having a characteristic progeni- 
tor mass set by fundamental physics, the la supernovae 
are not entirely homogeneous in their properties. Those 
found in star forming hosts are both more frequent (per 
unit stellar mass) and more luminous than their cousins 
in quiescent galaxies (Hamuy et al 1995; Scannapieco & 
Bildsten 2005; Sullivan et al 2006). 

In this Letter^ we look for similar trends linking the 
propertie of SGRBs and their hosts. We find an an- 
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ticorrelation between SGRB isotropic energy and host 
galaxy luminosity. We find no obvious selection effects 
that could produce such an effect, nor do we find evi- 
dence for any similar effect in a larger control sample of 
LGRBs. 

2. OBSERVATIONS 

We take our SGRB sample from the compilation by 
Berger (2008; hereafter B08) of aU SGRBs locahzed by 
the X-ray telescope (XRT) on the Swift satellite (Gehrels 
et al 2004). This sample contains 23 SGRBs, of which 12 
have reported redshifts. These 12 form the core sample 
for our study. They are further divided into 6 with op- 
tical afterglows, and hence the most accurate positions 
and most secure host galaxy identifications ("sample 1" 
of B08), and 6 with somewhat less accurate positions 
based X-ray data alone ("sample 2" of BOS). The two 
subsamples are similar in most properties (B08). 

We also form a comparison sample of 34 long bursts, 
by combining host galaxy information from Savaglio et 
al (2008) with burst durations from the GCN circulars 
and other GRB properties from Amati et al (2008). 

We calculate host galaxy absolute magnitudes as 
Babs = R- 51og(dL/10pc) + 2.51og(l + z). Here R 
is the observed R-band magnitude from B08, and di 
is the luminosity distance, calculated assuming Hq = 
71kms-iMpc"\ fio.m = 0.27, and f^o.A = 0.73. The 
resulting absolute magnitude corresponds to rest wave- 
length 6500A/(1 + z). For our SGRB sample, this al- 
ways between U and R band, and typically near B band 
(hence our notation). We calculate the isotropic equiva- 
lent gamma ray energy of the bursts directly from their 
fluence / as Eiso — 47rd|//(l + z). 

3. RESULTS 

The isotropic-equivalent energy of the SGRBs is corre- 
lated with the absolute magnitudes of their host galaxies, 
with a correlation coefRcient of 0.60. That is, the bright- 
est bursts tend to occupy the least luminous hosts, and 
vice versa (fig. [TJi). This is not expected from simple 
observational selection effects. If the samples are largely 
limited by flux and/or fluence, with most objects near the 
detection threshold, the range of distances in the sample 



should then lead to a positive correlation of burst energy 
with host luminosity. Such a positive correlation is seen 
for our control sample of long GRBs (figure [TJ)). 

The least luminous SGRBs are the closest, while the 
least luminous hosts span a wide range of redshifts. 
Gamma ray sensitivity limits inclusion in our sample 
more critically than does sensitivity to the host galaxy's 
starlight. If the host is not detected in the first relevant 
observation, more data can be taken, but if the GRB is 
not detected in the first observation, it can never enter 
the sample. 

If a low-luminosity event such as GRB 050509B oc- 
curred in a subluminous host at z < 0.5, we should still 
see it. We see perhaps one such event, GRB 070724, 
which has low isotropic energy and luminosity, but lives 
in a star forming host galaxy. We would also expect that 
if a highly luminous SGRB occurred in a large elliptical 
host galaxy, we should be able to observe such an event 
(and its luminous host) out to 2; > 0.8. We see no such 
events. This suggests that the more luminous SGRBs 
are associated with young progenitors of some type that 
is absent in old elliptical galaxies. 

We have tested the significance of the correlation in 
figure [T^ using a Spearman rank correlation test. The 
test yields a rank correlation coefficient of 0.50 and a 
significance level of i = 1.84 with 10 degrees of freedom. 
Formally, this corresponds to a 95% significance. 

However, the true significance is higher, since observa- 
tional selection effects work to produce the opposite sign 
of correlation. We have performed simple simulations 
of the correlation produced by observational selection in 
the absence of any true, physical relation between the 
burst and host properties. We first constructed indepen- 
dent probability distributions for -Eisoi Bats, and z. To 
these, we added selection probabilities based on gamma- 
ray fiuence and host apparent magnitude. Each of these 
probability distributions was specified by a simple, plau- 
sible functional form. The functions were adjusted until 
the simulations were able to reproduce separately the ID 
marginal distributions of /, Eiso, R, Babsi a-nd z seen in 
the data (all within the Poisson uncertainties). 

In each simulation, we constructed a "parent sample" 
of bursts, and applied the selection probabilities to sim- 
ulate a "selected sample," whose size is constrained to 
equal that of the observed sample. We then measure 
the correlation coefficient between Eiso and Babs in the 
simulation. We repeat this procedure many times, and 
compare the distribution of simulated correlation coeffi- 
cients to that measured in the observed sample. 

Among 10'' simulations, only 121 have coefficients 
above the observed value. Thus, by considering selec- 
tion effects, we estimate the significance of the SGRB 
energy - host luminosity anticorrelation at about 98.8%. 

While the simulations include several ad hoc functional 
forms, they capture the essential points of any realisitic 
model that does not have correlations between burst and 
host properties. First, the bursts and hosts span a wide 
range of intrinsic brightness. Second, the rate as a func- 
tion of rcdshift could be matched by the volume ele- 
ment modified by some modest amount of rate evolution. 
Third, the observational selection functions are mono- 
tonically decreasing over the range from the brightest 
to the faintest observed SGRBs and hosts. Fourth, the 
selection on the SGRB fiuence excludes many more sim- 



ulated events than does the selection on host luminosity. 
This induces a relatively weak correlation between host 
and burst brightness in the simulated samples. More 
strongly peaked intrinsic properties, or a more even bal- 
ance between rejection by either GRB or host properties, 
would tend to induce stronger correlations in the simu- 
lated sample. Thus, plausible changes in our simulations 
should not strongly affect our conclusions. 

Our results do depend substantially on GRB 050509b: 
If we exclude it from the sample, the remaining 11 
observed points yield a correlation coefficient of 0.28. 
This corresponds to about a 14% chance in the simu- 
lations. For comparison, a hypothetical sample of about 
40 SGRBs with a true correlation coefficient 0.28 would 
be significant at the same the 99% level as our full sam- 
ple. Provided that the identification of the GRB 050509b 
host is regarded as secure, the correlation does not de- 
pend critically on other bursts without optical transients. 
If we use the subset of SGRBs with optical transients 
(BOS's "sample 2") plus GRB 050509b (which has no 
OT but still a rather secure host detection) , the correla- 
tion coefficient becomes 0.63, and the simulations show 
a significance level of 95.6%. 

For reference, our control sample of LGRBs shows a 
correlation coefficient of -0.54 between isotropic energy 
and host absolute magnitude. We adapted our simula- 
tions to reproduce the marginal distributions of 2;, R, 
Babs, and so forth for the LGRB sample, and again com- 
pared the distribution of correlation coefficients to that 
observed for the LGRBs. Formally, this too yields a dif- 
ference at the 99% level between the simulations and the 
data, here in the sense that the data show a stronger pos- 
itive correlation between GRB and host brightness than 
do the simulations. Because the sign of the correlation 
matches the sign of the selection effects, and because the 
correlation is driven substantially by the X-ray flashes 
(020903 and 060218) at the bottom of the LGRB E^so 
distribution, we do not attach great physical significance 
to this result yet. However, it is intriguing and may merit 
further investigation in future. 

We also examined the relation between SGRB gamma 
ray luminosity and host galaxy star formation rate per 
unit stellar mass. Because the two most luminous SGRB 
hosts are early type galaxies, with tight upper limits on 
their emission-line derived SFR (Prochaska et al 2006, 
Berger 2008, Savaglio et al 2008), this plot clearly sepa- 
rates these two SGRBs from the rest of the sample. The 
specific star formation rate is defined as s = M^/M^,, 
where M^, is the star formation rate, and M^, the stellar 
mass of the galaxy. However, lacking the data to fit an 
accurate stellar mass to some galaxies, we use the quan- 
tity log(5) EE XogU^+QABabs = log (s)+l0g [M^/Lhost)- 

To demonstrate the utility of 5, we compared it to 
the published s (from Savaglio et al 2008) for the 34 
LGRB hosts from our control sample. We find log(5) sa 
0.41og(s X Gyr) - 7.9, with a scatter of 0.3 dex (RMS). 
The results (fig. [2^) show that the separation of the 
SGRB hosts into actively star forming galaxies and old 
stellar populations is accompanied by a separation of 
burst gamma-ray luminosity, with only the least lumi- 
nous bursts found in the old hosts. The corresponding 
plot for LGRBs (fig. ^Bp) shows no correlations beyond 
the sample selection requirement that the least luminous 
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Fig. 1. — The relation between isotropic equivalent gamma-ray 
energy and host galaxy absolute magnitude for short GRBs (first 
panel) and long GRBs (second panel). The observed correlation in 
the SORB sample is significant at about the 99% level (see text). 
The control sample of LGRBs shows a correlation consistent with 
that expected from observational selection. SGRB data points are 
labelled with the GRB ID number. Redshift ranges are denoted by 
color. For SGRBs, blue means z < 0.3, green means 0.3 < z < 0.7, 
and red means z > 0.7. For LGRBs, cyan means z < 0.3, blue 
means 0.3 < z < 1.0, green means 1.0 < z < 2.0, red means 
2.0 < z < 3.5, and magenta means z > 3.5. 



events be reasonably nearby. A comparison of the two 
"figures also shows the tendency (noted previously by, e.g., 
~B08) for even star-forming SGRB hosts to have less vig- 
orous star formation (i.e. lower S) than do LGRB hosts, 
_and also the tendency for SGRBs to have lower gamma 
-ray fluxes. 

4. DISCUSSION 

" While our sample is small, the anticorrelation between 
~SGRB Eiso and host luminosity is probably real, based 
on our Monte Carlo simulations. The similar anticorre- 
_lation between brightness of type la supernovae and the 
.stellar population ages of their host galaxies was first 
-described on the basis of a similarly small sample (13 su- 
-pernovae; Hamuy et al 1995). We now examine possible 
-interpretations for this result. 

Suppose, first, that two distinct mechanisms produce 
'SGRBs. Such a situation would arise naturally if two 
_or more of the different proposed progenitors actually 
_do produce SGRBs. Combinations that have been ex- 
-plored include coalescence of double neutron star bina- 
-ries and of NS-black hole binaries (Belczynski et al 2006; 
-O'Shaughnessy et al 2008; Troja et al 2008), and combi- 
nations of NS-NS binaries with giant flares from magne- 
tars (Chapman, Priddey, & Tanvir 2008). Additionally, 
a range of SGRB properties might arise from NS-NS bi- 
naries alone, provided those binaries evolved through a 
range of different binary interaction histories (Belczynski 
et al 2006; Salvaterra et al 2008). 

If one mechanism is associated with stars younger 
than 1 Gyr, and the other with ancient populations (age 
Hq ) , we can explain the properties of flgures [Iji and 
2^. The fraction of old-progenitor SGRBs occurring in 
quiescent galaxies should match the fraction of stellar 
mass in quiescent galaxies. The present, limited data 
shows two bursts in quiescent hosts, and a third simi- 
larly low-luminosity burst (070724) in a star-forming host 
galaxy. This matches well the fraction of stellar mass in 
elliptical galaxies at low redshift, which is 54% to 60% 
(Baldry et al 2004). Conversely, the fraction of more lu- 
minous SGRBs occurring in quiescent hosts is 0/9 in the 
present sample (fig. [2^). Under a two-population sce- 
nario, the mechanism responsible for the more luminous 
SGRBs should have a negligible rate for ages exceeding 
a few Gyr. 

If instead there is a single, continuous distribution of 
SGRB and host properties, it suggests that the SGRB 
progenitor's properties vary systematically with either 
age or metallicity, which vary systematically along the 
Hubble sequence and have direct effects on stellar-scale 
physics. 

If SGRBs are due to binary neutron star inspiral 
events, an anticorrelation of SGRB energy and the spe- 
cific star formation rate of the host galaxy could come 
from the dependence of inspiral time and available en- 
ergy on mass. Gravitational radiation gives an inspiral 
time of 
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(e.g.. Landau & Lifschitz 1958), where Mi and M2 are 
the masses of the two bodies, a is the semimajor axis 
of their orbit, and G and c are Newton's constant and 
the speed of light. Thus, if a star formation event yields 
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Fig. 2. — The relation between GRB luminosity and the specific 
star formation rate of the host galaxy, for short (first panel) and 
long (second panel) GRBs. The least luminous SGRBs are the 
two in massive early- type host galaxies. While some bursts in 
star forming hosts appear comparably faint (i.e., GRB 070724), 
we have yet to observe a bright burst in a quiescent host. For 
the long GRBs, there is no visible correlation of Lj and specific 
star formation rate. Point labels follow the same definitions as in 
figure 1. 



a population of compact object binaries with a range of 
mass, we expect a range of inspiral times. The SGRBs 
occurring in quiescent hosts would be those with low 
masses and/or wide initial separations. The available 
energy reservoir is i? ~ GMiM2/Rfinai, where R final is 
the characteristic size of the system when the GRB en- 
ergy is released. If the SGRB event produces a black 
hole, we expect Rfinai ~ Rgrav oc (Ml -I- M2), so 
that E ~ MiM2cV(Mi + M2). Then, for Mi « M2, 
tinspirai ~ QqE"^ . Characteristic stellar ages for SGRB 
hosts range from ^ 10 Gyr for quiescent hosts, to '^ 
s~^ ^ 1 Gyr for typical star forming SGRB hosts (see, 
e.g., Savaglio et al 2008). If this tenfold range of age cor- 
responds to a tenfold range of tinspirai 1 it would imply 
a range of ^ 2x in neutron star mass, which is compa- 
rable to the range of NS masses believed to exist in our 
Galaxy. 

However, the corresponding 2x range of energy cannot 
explain the range of Eiso observed in figure [T^. A larger 
range in SGRB energetics might follow if R final ~ Rns, 
especially for a neutron star equation of state that gives 
dRns/dM < 0. Moreover, uq may correlate with Mi and 
M2 in some complex way depending on on both binary 
star mass transfer and detailed supernova physics. While 
a more thorough understanding of the SGRB mechanism 
is clearly required before we can fully model the relation 
between host and SGRB properties, the existence of such 
a relation will provide valuable clues to the nature of the 
short GRBs. 

The correlation in figure [1^ can also provide a tool 
for estimating SGRB redshifts. \og{Eiso) and Babs are 
related empirically by log{Eiso/eig) ~ 50-|-0.725(i?abs + 
20.2). Defining V = (Il/VI + z, so that Eiso = (47r2?2)/ 
and Babs = Rab — 51og[(I?/10pc)], we can substitute for 
Eiso and Babs to obtain a distance estimate: 

-log(I?est/Gpc) - 0.012+0.129(i?Ai3-20)-0.1771og(/_6) 

(2) 
where /_6 is the SGRB fluence in units of 10 ^erg cm ^. 
The conversion from 'Dest to the estimated distance dL.est 
and redshift Zf,st must of course follow the cosmology 
we used to derive equation [2l For our primary sample 
of 12 SGRBs, the rms residual in \og{T)est) is 0.17 dex, 
and in redshift z{'Df.st) is 0.147. The residuals from the 
fit substantially exceed the uncertainties expected from 
fluence and host magnitude measurements, implying that 
most of this scatter is intrinsic. 

As an alternative, we can directly fit a relation of the 
form log(zest) = log(2o) + a{RAB - 20) + 61og(/-6) by 
minimizing /\zrm,s = {{zest - Zobs)'^ /{N - 1))^^^. This 
approach gives \og{zest) = —0.56 -I- Q.1Q7{Rab — 20) — 
0.141og(/_6), with Az™s = 0.142. This redshift esti- 
mator is easier to apply, and is formally independent of 
the adopted cosmological parameters, though it is less di- 
rectly linked to whatever physical mechanism links Eiso 
and Babs- Differences between the predicted redshifts 
from the two equations are all much smaller than the 
scatter in either. 

Another possibility is that the external medium around 
SGRBs plays some role in their luminosity. The inter- 
stellar medium is typically less dense in elliptical galax- 
ies than in actively star forming galaxies. However, this 
should only affect the observed flux from external shocks 
that GRB ejecta drive into the ambient medium. Ex- 



ternal shocks are believed to power the afterglow emis- 
sion, but not the GRB luminosity itself (e.g., Piran 2005). 
Since our observed anticorrelation involves the SGRB 7- 
ray flux itself, we disfavor this explanation. 

To conclude, we have found an anticorrelation between 
the energy of a short gamma ray burst and the luminos- 
ity of its host galaxy. Such an anticorrelation occurs with 
a probability of '^ 1% in simulations that account for ob- 
servational selection effects, and so is probably real. Its 
physical origin is unclear, though it is most likely due 
to a correlation between the age of an SGRB progenitor 
and the luminosity of the explosion. If this correlation 
is a continuous distribution, it can provide approximate 
redshift estimates for SGRBs. Time will tell whether this 



correlation holds for a larger sample, and whether it is a 
property of a single distribution, or if it instead reflects 
an underlying division of SGRBs into two physically dis- 
tinct sets. 
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